Frequency up-converted radiation from a cavity moving in vacuum 
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We calculate the photon emission of a high finesse cavity 
moving in vacuum. The cavity is treated as an open system. 
The field initially in the vacuum state accumulates a dephas- 
ing depending on the mirrors motion when bouncing back and 
forth inside the cavity. The dephasing is not linearized in our 
calculation, so that qualitatively new effects like pulse shaping 
in the time domain and frequency up-conversion in the spec- 
trum are found. Furthermore we predict the existence of a 
threshold above which the system should show self-sustained 
oscillations. 

PACS: 42.50.Lc - 03.70. -Hk - 12.20.-m 

Vacuum field fluctuations exert radiation pressure on 
scatterers in vacuum. For a pair of mirrors at rest this 
effect is well known as Casimir effect When a mir- 
ror is moving radiation pressure of vacuum fluctuations 
leads to a dissipative force which opposes itself to the 
mirrors motion. This force is known to arise as soon as 
the mirror has a non- uniform acceleration Accord- 
ingly the electromagnetic field remains not in the vac- 
uum state but photons are emitted by the scatterer into 
vacuum . Radiation from a moving mirror and the as- 
sociated radiation reaction force imply that dissipative 
effects are associated with the motion of mirrors in vac- 
uum, although this motion has no further reference than 
vacuum itself. Since these effects challenge the principle 
of relativity of motion in vacuum, it would be very im- 
portant to obtain experimental evidence for them and to 
study their characteristics in detail. 

Motion-induced radiation can be interpreted as a result 
of dephasing of vacuum fields depending on the mirrors 
motion. The order of magnitude of the dephasing is ex- 
pected to be the ratio between the mirror's velocity v and 
the speed of light c. For most conceivable motion of a 
macroscopic object, the velocity v cannot greatly exceed 
the sound velocity and is thereby much slower than that 
of light. This is why motion- induced radiation is very 
small for a single mirror oscillating in vacuum. This con- 
clusion holds for perfectly reflecting mirrors as well as for 
partly transmitting ones. 

A number of works have been devoted to photon pro- 
duction inside the cavity built with a pair of perfectly 
reflecting mirrors moving in vacuum H-R]. However no 
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predictions can be made for the amount of radiation emit- 
ted outside the cavity when the resonator is treated as 
a closed system. In contrast the resonant enhancement 
is found to be determined by the cavity finesse when the 
cavity is treated as an open system from which the 
photons can escape. Motion induced radiation, that is 
photon emission outside a cavity oscillating in vacuum 
is resonantly enhanced by the cavity finesse when com- 
pared to the radiation from a single oscillating mirror. 
The resonant enhancement occurs when the mechanical 
frequency is a multiple of the lowest cavity mode. Even 
and odd multiples correspond respectively to breathing 
modes, where the mechanical cavity length changes pe- 
riodically, and to translation modes, where the cavity 
moves as a whole Q . The latter effect reminds radiation 
from a single mirror inasmuch as vacuum fiuctuations are 
the only reference for the cavity motion. However the or- 
der of magnitude of photon emission may greatly exceed 
the one from a single mirror. From an experimental point 
of view the cavity is so far the most interesting system to 
look for an experimental observation of dissipative effects 
of vacuum fiuctuations. 

Inside a cavity the field undergoes many reflections be- 
fore leaving the cavity through one of the mirrors. The 
number of round-trips of the field is roughly given by the 
cavity finesse. In loose terms, one may define an effective 
velocity where the physical velocity v normalized by the 
speed of light is multiplied by the number of round-trips 
inside the cavity. Effective velocity and thus motion- 
induced radiation become the larger the higher is the 
cavity finesse. The effective velocity is no longer a ma- 
terial velocity and may therefore approach the speed of 
light. In contrast to the single mirror's case, qualitatively 
new effects are expected, such as the formation of a pulse 
bouncing back and forth in the cavity Since the 

pulse duration is shorter than the time of fiight through 
the cavity, the radiation spectrum should also contain 
various frequencies corresponding to higher-order cavity 
modes and thus exceeding the mechanical frequency. 

These effects cannot be obtained with a linear treat- 
ment as the one used in |^ . In such a treatment the field 
scattering is supposed to be linear in the mirrors motion. 
The field-mirror interaction corresponds to a coupling of 
the vacuum field radiation pressure, which is quadratic in 
the field, to the mirrors mechanical motion. Photons are 
thus emitted in pairs. In the linear regime the generation 
of motion-induced radiation is analogous to a parametric 
process during which the mechanical excitation is trans- 
formed into a pair of photons. Due to energy conserva- 
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tion the sum of their frequencies equals the oscillation 
frequency. Therefore motion-induced photons are only 
emitted at frequencies smaller than the excitation fre- 
quency. The linear assumption is correct as long as the 
total field dephasing due to interaction with the moving 
reflector remains small. The field dephasing from one re- 
flection scales with the mirrors velocity over the speed 
of light. The linear assumption is always satisfied for a 
single macroscopic mirror. 

However for a cavity the crucial parameter is the ef- 
fective velocity and the total field dephasing can become 
important for a high finesse cavity. As a consequence 
we expect frequency multiplication to occur which gen- 
erates frequencies larger than the mechanical excitation 
frequency. The linear treatment, which predicts the emis- 
sion frequencies to be smaller than the oscillation fre- 
quency, then looses its validity and has to be replaced by 
a treatment which fully accounts for a large field dephas- 
ing produced through successive reflections of the field 
onto the mirrors. This treatment will be called non-linear 
in the following although the scattering is still linear in 
the field. 

The aim of the present paper is to give a treatment 
of the radiation emitted by a cavity moving in vacuum 
which takes into account both the effect of accumulated 
field dephasing and the open character of the cavity. In 
particular we will study the explicit dependence of mo- 
tional radiation on the two experimentally accessible pa- 
rameters, the mirrors velocity and the cavity finesse. We 
will first introduce general calculation techniques for a 
single mirror and for a cavity moving in vacuum with 
an arbitrary motion. We will then concentrate on a par- 
ticular motion, the harmonic oscillation of the mirrors, 
which allows to evaluate in closed analytical form the 
correlation functions of the radiation through a special 
parametrization of the motion. We will give expressions 
for the time-dependent radiated energy as well as the 
frequency-dependent radiation spectrum. 



I. SINGLE MIRROR MOVING IN VACUUM 

Neglecting all effects related to polarizations, the elec- 
tromagnetic field $ is considered as a scalar function of 
one space variable, t and x, and will be written as a sum 
of two counterpropagating components y and ip which 
are functions of two light-cone variables u and v 



U — t — X V = t + X 



(1) 



For the sake of simplicity, the velocity of light is set to 
unity. We limit ourselves here to two-dimensional space- 
time calculations. As is well known from the analysis 
of squeezing experiments px|], the transverse structure 
of the cavity modes does not change appreciably the re- 
sults obtained from this simplified model. Each trans- 
verse mode is correctly described by a two-dimensional 



model as soon as the size of the mirrors is larger than the 
spot size associated with the mode. The two-dimensional 
model thus corresponds to a situation where one trans- 
verse mode is efficiently coupled to the moving mirrors. 

We now represent the mirror's motion cither by a func- 
tion X — q(t) associated with the trajectory or as a 
monotonous function v = V{u) relating the light cone 
variables u and v of the light rays intersecting on the 
mirror's trajectory. A propagation component tpout of 
the output field can then be written as a function of the 
input fields (pin and tpin and of the function V 



V'out = VRljjin oV + Vrifiir 

T =1- R 
{hog) (u) = h{g{u)) 



(2) 



The symbol o represents the composition law for func- 
tions. The coefficients VR and Vt are the reflection 
and transmission amplitudes describing scattering upon 
the mirror. For the sake of simplicity, we have assumed 
these coefficients to be real and frequency independent. 

We now recall the standard calculation of the energy 
radiated into vacuum by the moving mirror ||^. The in- 
put fields are supposed to be in the vacuum state and 
characterized by the correlation function 
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In |u — u| — Tr^iu ~ u) (3) 



The first term corresponds to the anticommutator and is 
state-dependent whereas the second term gives the com- 
mutator contribution and is state- independent. Since e 
is the sign function, it is clear that the field commutator 
remains unchanged under the transformation (^) where u 
is replaced by a monotonous function V {u). The change 
of the correlation function between the input and output 
field is given by the following function which depends on 
the field anticommutators only 



fiR I V'{u)V'{u) 1 
~ ~ 4^ \ {V{u) - V{u)Y ~ {u - uf 



(4) 



Throughout the paper, the prime signifies a derivative 
of a function with respect to its argument. The energy 
density e„(u) radiated per unit time is given by the func- 
tion C{u, u = u) evaluated at coinciding points through 
a point splitting regularization procedure S 



TiR 

eu{u) = C(u, u) = ~——SV {u) 



sv 



V'" 



v 



(5) 



The function SV is the Schwarzian derivative of V . No 
radiation is emitted when the reflector has a uniform ac- 
celeration, which corresponds to a vanishing Schwarzian 
derivative SV . The total energy radiated by the moving 
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mirror can then be obtained by integrating the energy 
density over u. In the following, we will concentrate on 
the particular case of an oscillatory motion of mechani- 
cal frequency O. In this case the energy E^^ radiated per 
period is read as 

IT 

Eu^ eu{u)du (6) 
Jo 

In order to characterize the radiation we have also to 
describe its spectral properties. The radiation spectrum 
may be represented as a density of photons obtained from 
the Fourier transform of the two point function C{u, u). 
We will turn to its description later on. 

II. CAVITY MOVING IN VACUUM 




X 

FIG. 1. Space-time diagram of on arbitrary field trajectory 
bouncing back and forth inside the cavity. Both mirrors are 
supposed to follow a harmonic motion. Light rays are indi- 
cated by null lines, i.e. straight lines making a 45°-angle with 
the time and space axis. 



The vacuum field is defined on both sides of the cavity 
as in the previous section. The relation between input 
and output fields is a generalization of (j^) which corre- 
sponds to the standard Fabry-Pcrot theory 

"^out = - V^V'in oV-i + \fR\T2 ^ r"-0in O 14 

ri>0 

ji>0 

Ti = 1 - i?i T2 = 1 - i?2 

r = /rT/R^ = e-2p (7) 

The reflection and transmission amplitudes of the two 
mirrors are related through unitarity conditions. The co- 
efficient T determines the attenuation factor of the field 
on a single cavity round-trip. It can also be written as a 
function of the cavity losses p. Throughout the paper we 
will use p when we consider the experimentally interest- 
ing case of a high finesse cavity with p In the more 
general case the reflection coefficient r will be used. The 
functions C/„ and Vn represent the light cone variables 
associated with the various input rays which are trans- 
formed into the output light ray u by the cavity. When 
the cavity is at rest they are given by simple relations 

Vn (u) = u — 2nL 

Vn (u) =u- {2n + 1)L (8) 

The length L of the cavity is measured as a time of flight 
between the two mirrors; the two mirrors are supposed 
to be located at x = ±L/2 respectively; the ray V^i 
represents the particular case where the field has been 
directly reflected back by the flrst encountered mirror 
without entering the cavity. 

We may deduce the held emitted by the vibrating cav- 
ity through the same expression (|7|) as for the motionless 
cavity, but with functions Un and Vn now given by the 
procedure sketched on flgure 0. 



The functions [/„ and Vn corresponding to the various 
light rays in figure |l| are built up through a functional 
iteration, 

Un — f2n Vn = /2n+l 

/_! - .9 /o - / 

f2n = 9^^ ° /2n-l /2n+l = h O f2n (9) 

where / is the identity function {I (u) = u), while the 
two functions h and g"^ (inverse of g) represent the tra- 
jectories of the mirrors 

X — qi{t) =^ V = h{u) 

X ^ q2{t) =^ V = g{u) ^ u = g^^ (v) (10) 

The function f2n results from n successive compositions 
of the function g~^ o h. This construction is quite anal- 
ogous to the one described in Q, but it is written here 
such that it may be applied to an open cavity. 

In a linear treatment |^ the total field dephasing re- 
mains small. At every reflection the fleld acquires a de- 
phasing due to the mirrors motion 2qi{t), (i = 1,2). In 
this case the composition of functions (g) is reduced to 
the summation of the mirrors motion. The accumulated 
dephasing after M roundtrips inside the cavity is now 
seen to be simply Af times the dephasing due to a single 
roundtrip. The scattered field then has a temporal vari- 
ation which reproduces the mirrors motion. With these 
approximations the results of the linear treatment are 
recovered. 

However, if the fleld undergoes a great number of 
roundtrips inside the cavity, the total dephasing does 
not remain small and a linearization is not valid any- 
more. The dephasing has then to be calculated through 
the general composition law (^ . Following the same pro- 
cedure as for the single mirror, we calculate the density 
of energy e„(M) radiated to the right through the two- 
point correlation function deflned as in (^) by letting u 
come to coincidence with u. As all functions now depend 
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on a single parameter u, we omit this parameter in the 
expression for the energy density 
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I 6 (/-I - /2„+l) 
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n>0 



n>0 



T1T2 ^ r" 



/2n/2 



J2nJ2 



{fin ~ /2m)^ 

» fln+lflm+l 
(/2ri+l — /2m+l)^ 



(11) 



Compared to the radiated energy density of a single mir- 
ror we now find a sum of Schwarzian derivatives cor- 
responding to different numbers n of round-trips inside 
the cavity, as well as new terms arising from the inter- 
ference between light rays having undergone a different 
number of roundtrips. 

The derivatives appearing in the upper equation are it- 
eratively deduced from each other through relations (||) 
and the chain rules associated with derivation of com- 
posed functions 



{g o h)' = h' g' o h 
S {g oh) = Sh + h'^ (Sg) o h 



(12) 



In the general case of arbitrary motions of the two mir- 
rors, the various relations which have been written in the 
present section allow to compute the energy density radi- 
ated by a cavity built with partly transmitting mirrors. 



III. HARMONIC MOTIONS AND PERIODIC 
ORBITS 

From now on we focus our attention on configurations 
which have been shown to be the most efficient ones to 
generate motion-induced radiation ||^ and which further- 
more allow us to put the problem in a simpler form. 

We consider that the two mirrors follow harmonic mo- 
tions at such a frequency that the motion-induced effects, 
i.e. motional radiation and motional force, are resonantly 
enhanced by the multiple interference taking place inside 
the cavity. The frequency f2 of the harmonic motion is 
thus supposed to be such that QL is a multiple of tt 



flL ^ Kn 



(13) 



The amplitudes of the two motions are supposed to have 
the same absolute value with either opposite or identical 
signs depending on the parity of the integer number K 



nqi{t) 



Kn 



/3 sin nt 



{K + l)7r 



nq2{t) 



Kit 
/3 = th(Q;) 



/3sin VLt 



(i^ + l)7r' 



(14) 



We have written the two equations of motion in terms of 
dimensionless numbers. In particular, one distinguishes 
two parameters a and /3. /3 represents the ratio between 
the maximal velocity of the mirrors and the velocity of 
light while a plays the role of a rapidity which will be 
found to add up through successive reflections when one 
considers the composed motion of both mirrors. 

When K is an even number, the upper equations de- 
scribe a situation where the two mirrors are oscillating 
such that the length of the cavity changes periodically. 
In the opposite case, the cavity performs a global oscil- 
lation with its length kept constant. The two cases will 
be called even and odd modes respectively in the fol- 
lowing. A cavity motion corresponding to an odd mode 
is reminiscent of motion-induced radiation from a single 
oscillating mirror, but in addition here the radiation is 
resonantly enhanced inside the cavity. These statements 
follow from the linearized approach but we expect a 
similar behavior to take place in the full treatment de- 
veloped in the present paper. 

There exist periodic orbits such that the optical length 
seen by the field bouncing back and forth in the cavity 
is the same on successive round-trips despite the motion 
of the mirrors . These orbits correspond to particular 
values u of the light-cone variable u such that the iter- 
ation procedure leads to expressions similar to the ones 
obtained when the mirrors are at rest (cf. (P). Although 
definition (^ of /p is different depending on whether p 
is even or odd, the periodic orbits generalize the usual 
resonance condition of the Fabry-Perot theory 



(15) 



They therefore give rise to a constructive interference ef- 
fect, analogous to that occurring for a motionless cav- 
ity. Since the light rays corresponding to a periodic or- 
bit encounter the mirrors at the same position after an 
arbitrary number of round-trips, the composition (|^) of 
motions leads to a simple power law for the derivatives 
evaluated after n roundtrips as well as for the Schwarzian 
derivatives 



5/p(u)=5/({t)y 



(16) 



There exist two sets of periodic orbits which correspond 
to opposite values of a. in ([l6|). These two periodic or- 
bits attract (respectively repel) the neighboring trajec- 
tories, when e*" is greater (respectively smaller) than 
unity. Only the attractive orbit is expected to give rise 
to a large enhancement of the motional radiation. 

In expression (^) of the energy density the Schwarzian 
derivative (^) is multiplied by the squared reflection 
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coefRcient after n roundtrips r^". Summation over the 
number of roundtrips leads then to a geometric progres- 
sion of 7-^"e^"". The first factor represents the atten- 
uation of the energy density associated with the cavity 
losses through the two mirrors. The second one accounts 
for the parametric amplification of the field associated 
with the mirrors' motion. As a consequence, the energy 
density takes large values when the parametric amplifica- 
tion compensates the losses. In fact, a divergence of the 
energy density should occur when re''" approaches unity. 
This corresponds exactly to the oscillation threshold of a 
mechanically excited parametric amplifier. Let us notice 
that the approach developed in the present paper does 
not remain valid above this threshold. 

We have focussed our attention here on the case where 
periodic orbits correspond to light rays meeting the two 
mirrors at their mean positions, respectively — and -j. 
There exist more general situations where the light rays 
meet the mirrors at other positions which will not be 
considered here. Notice that the particular case studied 
in the present paper is interesting from an experimental 
point of view since it corresponds to the maximum value 
of the parameter a for motions having a given frequency 
and a given amplitude. 

In the following, we will restrict our attention to the 
cases of practical interest where the physical velocity of 
the mirror remains small when compared to the veloc- 
ity of light. As discussed in the Introduction, this con- 
dition is always met for macroscopic mirrors. It implies 
that a single refiection produces a small dephasing on the 
field and, thereby, small radiation effects. Precisely, this 
means that the quantity Sf which appears in ( |l6|) has an 
extremely small value while the factor e^" is very close 
to unity. It follows that a large number n of roundtrips is 
needed to obtain a factor e^"" differing appreciably from 
unity and therefore giving rise to a noticeable radiation. 
We show now that this assumption permits to perform 
the functional iteration in an analytical manner. 

The crucial point is that the functional iteration (|p) 
may in this case be restricted to the sub-space of peri- 
odic functions h corresponding to homographic relations 
between the phases e*^" and e*^^''^") 



ae 



m 



(17) 



a and b are two complex constants and a* and b* their 
complex conjugates which can be gathered in a matrix 
Aih) associated with the function h. Attention may 
be restricted to matrices having a determinant equal to 
unity. In the sub-space of functions (p7|), the composi- 
tion of two functions merely corresponds to a product of 
their corresponding matrices. 



A{hog)^A{h)A{g) 



(18) 



Rigorously speaking, the function ( p^ ) does not corre- 
spond to sinusoidal trajectories ( |l^ but rather to specific 



trajectories already considered by Law for dealing with 
photon production inside a closed cavity 



sin (f2(7 — (/)q) — —/3 sin (rJi — (/)(,) 



(19) 



where the reduced velocity (3 gives the mirrors velocity 
compared to the speed of light. Variations of the phase 
factors 0a and (j)b amount to displacements of the trajec- 
tory along the space and time axis. At the limit of small 
velocities however, the trajectory (|l^) is reduced to an 
ordinary sinusoidal motion (n4h 



(20) 



13 <^ I =^ Vtq = 4>a~ 13 sin {fit - 4>b) 



The difference between the two motions scales as the cube 
of velocity and is therefore extremely small for real- 
istic motions of macroscopic mirrors. The effect of the 
trajectory (|l^) is thus indistinguishable from the effect 
of the sinusoidal motion for a single refiection. 

If one considers two oscillating mirrors, the two func- 
tions h and g, corresponding to the first and second mir- 
ror, are associated with two matrices A{h) and A{g) 
respectively. The matrix components are chosen to fit 
equations (Ol) of motion of the two mirrors 



A{h) 



{—i)^cha z^+'sha 
(— i)-^+-'^shQ; i^cha 

i^cha (— i)^+'sha 



(21) 



As the composition of the two functions corresponds to 
a product of their matrices, the composition law natu- 
rally produces a homographic function when the number 
of refiections becomes large. This essential feature will 
be used in the following to compute the temporal and 
spectral characteristics of motional radiation from purely 
algebraic manipulations of the associated matrices. In 
particular, the functions fp in equation ^ correspond- 
ing to successive refiections of the field inside the cavity 
are obtained through matrix multiplication 



A2K+1 



{-i)^Pshpa \ 



(22) 



The discrepancy between the composed functions built 
on the two motions (14) and ( [l9| ) does not affect the 
results if /? ^ 1. More precisely, the difference between 
the composed functions fp built on the two motions ([l^ ) 
and (|l^) remains of the order of when the number of 
iterations increases. 



IV. PULSE SHAPING 

In the following, we will analyze the case of a single 
mirror following a trajectory ( p^ ) with an arbitrary ve- 
locity parameter a. Although the hypothesis of a large 
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rapidity a is not realistic for a single mirror, it can be 
used as a model for the composition of a large number of 
round trips inside the cavity. We will then come to the 
full treatment of the cavity where interferences have also 
to be accounted for. 

The derivative h'{u) may be written from ( |l7|) as 



h'{u) 



1 



he 



'iQ,u I 



(23) 



h'{u) oscillates between the extremal values exp (±2a) 
which correspond to physical velocities of the mirror ±/3. 
For homographic functions ([l^), the Schwarzian deriva- 
tive has the simple form 



Sh=-il 



(24) 



The total energy radiated to the right by a single mirror 
is then obtained by averaging the energy density (|^) over 
one oscillation period 



nm 2 



(25) 



This energy does not saturate when the parameter a in- 
creases although the velocity scales as tha and remains 
smaller than the velocity of light. The radiated energy is 
always greater than the squared velocity th^a which was 
the value suggested by the linear treatment ||] . However 
it is impossible to obtain appreciable radiation with a 
single oscillating mirror for velocities small compared to 
the speed of light. In the realistic case of a mirror mov- 
ing at a small velocity the radiated energy as well as the 
number of emitted photons scale with . 

We come now to the energy density radiated by an os- 
cillating cavity, assuming that a large number of round- 
trips is necessary to compensate the small velocity of the 
mirrors and thus get a noticeable radiation. The energy 
density e„ may be obtained from (11) by using the fol- 
lowing properties of the Schwarzian derivative of fp and 
of the first-order derivative fp respectively 



Sfp 

f 
jp 



ch2pa + {—l)^sh2pa sin flu 



(26) 



To plot the energy density for different linear and non- 
linear regimes we introduce effective quantities 



Qeff = 2a/p 
/3off = th(acff) 



(27) 



The effective rapidity aos is given by the roundtrip value 
of the rapidity multiplied by the cavity finesse In 
contrast to the mechanical velocity v normalized by the 
speed of light, which has to remain much smaller than 
1 in any physical situation, the corresponding effective 



velocity f3eS can become an important fraction of the 
speed of light when the field undergoes a large number 
of reflections inside the cavity. The maximal value of 
the effective reduced velocity is limited to /3cS ~ 0.76 
by the divergence of the energy density at acs — 1. This 
value corresponds indeed to the threshold re'*" = 1 which 
has already been mentioned previously for the periodic 
orbits. 

The variation of the energy density for different param- 
eters Qfeff is presented on figure ^. In the linear regime 
where aeff ^ 1 the temporal variation of the emitted en- 
ergy is sinusoidal. When a^s increases the energy concen- 
trates in pulses which are periodically emitted by the cav- 
ity. This pulse shaping becomes the more pronounced, 
the width of the pulses the smaller, the larger becomes 
the effective rapidity and thus the accumulated field de- 
phasing. 

In the same manner as for the single mirror the total 
energy radiated by the cavity is computed by averaging 
the energy density over one oscillation period 27r/f2. 




FIG. 2. Energy density emitted to the right by the cavity 
as a function of time for different effective rapidities Oeff and 
r — 0.99. The top line of the frame corresponds to 10~^hQ^ . 
With increasing values of the effective rapidity the energy 
starts to concentrate in pulses emitted periodically by the 
cavity. The pulses become very sharp close to the threshold 
of oscillation Qeff ~ 1. 

As previously we will restrict our attention to a realistic 
case, where the mirrors velocity is small which justifies 
the use of the homographic relations. In this case we 
may replace the differences of functions in the denomina- 
tors of (11) by their motionless values (|§). Higher order 



contributions decrease with the inverse of their squared 
order. The sums over the number of roundtrips can then 
be performed and we find the following expression for the 
energy radiated to the right 



nnR2 ,2 , nnT2 

£/„ = sh a H 
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(Cn(a) + C«(-a)-2) 
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nnT2 



E : 



C«(a) 



(1 -e-4'')e2«+Ti(l 



o4(q-p) 



1=1 



(28) 



The energy radiated to the left can be obtained from 
the above formula by interchanging the indices 1 and 2 
of the reflection and transmission coefficients. The total 
energy radiated in both directions is then evaluated as 
the sum of the two contributions 



E — E,, 



Ey 



(29) 



The total radiated energy would diverge for a = p. How- 
ever this limit is not reached as the energy density al- 
ready diverges when ia/ p approaches unity. The reason 
for this difference is simply that the widths of the var- 
ious contributions to the energy density decrease when 
the number of roundtrips n increases, so that the con- 
tribution to the integrated energy increases less rapidly 
than the peak value of the energy density. As a conse- 
quence, the divergence of the peak density occurs before 
the divergence of the integrated energy. We remind here 
that our calculation does not remain valid above the os- 
cillation threshold. Still it may be expected that a large 
amount of radiation is emitted above threshold. 

For experimental reasons one might also be interested 
in the amount of energy stored inside the cavity. Let us 
first remark that inside the cavity the fields propagat- 
ing to the right and to the left are not independent from 
each other due to the boundary conditions. All intracav- 
ity quantities thus contain the field components of both 
directions and are obtained as the sum of the two contri- 
butions. Having this argument in mind the intracavity 
energy, integrated over the cavity length L — Kn/Vl, is 
then with the same notations as above found to be 

^ = ^(C(") + C(-«)-2) 



- (^(a)C(«) + ^(-«)C(-") - 2^(0)) 

C(a) = ^(C«(a) + C.(a)) 



(30) 



The energy density is here directly expressed with respect 
to the static Casimir energy which is recovered when the 
cavity is motionless . This result is due to the fact that 
the vacuum outside and inside the cavity is not the same 
but differs exactly by this amount of energy. 

In order to obtain an appreciable value for the radiated 
energy if the cavity is moving at a small velocity, it is nec- 
essary to consider a high finesse cavity p <C 1 , keeping in 
mind that the finesse should be limited by the condition 
aeff < 1- Using these assumptions equations (|2^j3^) may 
be approximated as follows by expanding separately the 
common denominator and numerators 



24 
P 



K- 



1 



a 



1 

2 



pa 



Kjp^ 



a<'-«l 



(31) 



The first term in the radiated energy is due to the field 
which is directly reflected by the two mirrors without en- 
tering the cavity. This term corresponds to the expres- 
sion for motion-induced radiation from a single perfectly 
reflecting mirror. The second term has its origin in the 
field which has traversed the cavity and thus accumu- 
lated a much more important dephasing than the singly 
reflected field. Neglecting in the denominator leads to 
the linear result presented in The linear approxima- 
tion is found to be rigorously valid if the rapidity is much 
smaller than the cavity losses (a <C p). Furthermore the 
present treatment allows us to calculate motion-induced 
radiation when the fleld dephasing becomes large due to 
accumulation on a large number of reflections. Equa- 
tions (|l|) have a range of validity extending up to the 
threshold acff = 1- 

The intracavity energy contains only the term corre- 
sponding to the field which has entered the cavity. Its 
expression can also be deduced from the radiated energy 
through a detailed balance argument, which goes as fol- 
lows The energy inside the cavity can be obtained 
from the radiated energy by multiplying it by two fac- 
tors f2/(27r) and 2L/{Ap). The first factor is due to the 
fact that the radiated energy is the energy density inte- 
grated over one oscillation period. During one roundtrip 
of duration 2L each photon has the probability Ap - cor- 
responding to the energy transmission coefficient of the 
two mirrors - to escape from the cavity. The present non- 
linear evaluation of the intracavity energy gives indeed 
the same result as the balance argument in the limiting 
case of a high finesse cavity (cf. equations (^l|)). How- 
ever, as equations (|28| , |30| ) show the balance argument is 
not true for a cavity with arbitrary refiection and trans- 
mission coefficients. A remarkable consequence of equa- 
tions (^Tj) is that the non-linear calculation is necessary 
as soon as the number of photons inside the cavity be- 
comes of the order of unity. 

Interesting remarks can be made concerning the partic- 
ular case K = 1. Clearly equations ( pl| ) show that for a 
high finesse cavity no enhancement of photon production 
inside the cavity can be obtained when the mechanical ex- 
citation frequency equals the lowest cavity mode {K = 1) 
in accordance with results in ref. H,^. In this case the 
energy E radiated by the cavity corresponds to the one 
emitted by a single mirror and the motional intracavity 
energy vanishes. However, in the general case of arbitrary 
cavity finesse motion-induced photons are also found for 
the K ~ 1 mode (cf. equation (pO|)). The key point is 
that as far as classical light rays are concerned the mode 
K = 1 behaves like all other modes ||^. However the 
field dephasing and thus motion-induced radiation is not 
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only determined by the behavior of the hght rays but also 
by the cavity which plays the part of a filtering function 
and suppresses photons at zero frequency. As a conse- 
quence of the coupling to radiation pressure photons are 
not emitted singly but in pairs. Thus motion-induced 
radiation is enhanced by the cavity if all photons are 
emitted into a cavity mode, the sum of their frequencies 
being equal to the mechanical oscillation frequency. In 
order to fulfill this condition when the cavity oscillates 
with the frequency of the lowest cavity mode photons 
have to be emitted at zero frequency. The cavity sup- 
presses those photons the more efficiently the higher is 
the cavity finesse. Thus motion-induced photons for the 
K — 1 mode can be found in the bad cavity limit but not 
in the high finesse limit. 

Coming back to the general case, we emphasize that 
equations (^ij) remain valid up to the threshold p = 2a, 
when the cavity finesse is increased with the ampli- 
tude of motion kept constant. Below this value motion- 
induced radiation is amplified inside the cavity, but the 
cavity losses exceed the amplification gain. As expres- 
sions (^l]) are monotonous in p their maximum values are 
thus reached at threshold. In this regime we then find a 
maximum of the radiated energy which depends linearly 
on p. Comparing this value to the maximum energy emit- 
ted by a single oscillating mirror shows a gain of the order 
of the cavity finesse by considering a cavity instead of a 
single mirror. The cavity is thus a much more favorable 
system to produce motion-induced radiation than a sin- 
gle mirror. Furthermore if one increases the cavity finesse 
above its threshold value the roundtrip amplification of 
the field due to the mirrors motion should exceed the 
cavity losses and the system should enter a regime of ex- 
ponential amplification. Without further calculations we 
then expect the oscillating cavity to emit photon pulses 
of much higher intensity above threshold than below. 



V. FREQUENCY UP-CONVERSION 

We now turn to the calculation of the radiation spec- 
trum where we will proceed as previously by first study- 
ing the case of a single moving mirror and afterwards the 
one of the oscillating cavity. 

The scattering field equation (0) writes in Fourier space 



„2iPQ(«) 



n 



to 

n 



(32) 



where we have introduced the reduced frequencies v and 
v normalized with respect to the mechanical frequency 
fi. The field dephasing of the output field is determined 
by V{u) and thus associated with the mirrors position 
Q{u) which is easily calculated from (|l7| ) 



nq {u) = n 



Viu) 



arctg 



f3 cos flu 
1 + P sin Qu 



(33) 



is a periodic function and can thus be developed 



into a Fourier series with discrete coefficients 7m 

e-^'''^Q(^) = 7™[i^]e-'""^^" (34) 

m 

which will determine the radiation spectrum. If the 
mirrors motion were sinusoidal these Fourier coefficients 
would be given by Bessel functions of different orders. We 
have now to evaluate these coefficients for a homographic 
trajectory (19). 

To this aim, we first write the field dephasing (|3^ ) as 



I + if3e 



-inu \ 



1 - i/Se'^" 



(35) 



The Fourier coefficients (p9) may be rewritten in terms 
of an hypergeometric seriesn 



7m H = (-i)"^^ - sin (vrzy) G„ (i', /?) 



l>0 



r(m-f-i + o 



n 



(36) 



The radiation spectrum, that is the spectral density of 
the photon number per unit time and defined for positive 
frequencies, is then given by M 



riu = R 



-\l-m\m - v] 



7n>l> 
.2/ 



i?^^^ V.(m-.)|G„.(.,/5)r 



(37) 



The total energy may be recovered as the integral of the 
spectrum as well as the integral of the energy density. 
The radiation spectrum vanishes for all values v equal 
to a natural number, that is for all frequencies lo equal 
to a multiple of Q. The spectrum thus decomposes into 
a succession of arches, each limited by two successive 
multiples of the excitation frequency. 

The upper expression corresponds to reflection upon a 
single moving mirror. In order to get the spectral distri- 
bution of radiation from the cavity, we now have to take 
into account the interferences between light rays having 
undergone different number of reflections inside the cav- 
ity. We proceed in the same manner as in the case of 
a single moving mirror by splitting the function fp into 
two parts. The first part is linear in the parameter u 
while the second one Qp is induced by the motion and 
harmonic like the mirrors' motion 



^The coefficient Gm{i^,P) is directly related to the hyperge- 
ometric function F(i^, m — u;m + 1; P^) defined for instance 
in We have also used property 8.334 of the F-function in 
the same reference. 
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fp (u) = u-pL + 2Qp (u) 
ilQp[u) = arctg ' 



1 + /3p sin flu 



/3p = (-l)^th(pa) 



(38) 



The round-trip dephasing 2L corresponds to the case of 
periodic orbits so that the function Qp vanishes at these 
points. The periodic function e^^'^^p can be developed 
into a Fourier series with cocfhcients now depending on 
the number of round trips p. We proceed as previously 
to find the Fourier coefficients 



27r 



+2, -Sin iTTiy)iKnuv 



G„,{v,l3p) (39) 



and the radiation spectrum 
sin^ {ttv) 



/ TD rji \ n —2i7rKi^(n-\-l)/^ (a \ 



n>0 



(40) 



Let us mention that we recover the predictions of the lin- 
earized treatment for a motion with a small velocity by 
keeping only the lowest-order term m = 1 in the hyperge- 
ometric series. The spectrum is then parabolic and found 
to be restricted to the frequency range corresponding to 
the first arch ||]. 

Figure |^ shows the radiation spectrum for an effective 
rapidity a^s = 0.9 near the threshold of parametric oscil- 
lation. The spectrum shown here is plotted for a cavity 
oscillating globally at a frequency of fl = 3tt/L. This 
means that the cavity performs three oscillations during 
one roundtrip of the field inside the cavity. The dashed 
line was obtained by putting formally K = in equa- 
tion (^0|) which eliminates the phase factors responsible 
for the interferences. It may be interpreted as the spec- 
trum of radiation emitted by a single oscillating mirror 
averaged over the effective velocity. Clearly photons can 
be created by higher-order harmonics of the motion as 
well as by the fundamental one as soon as the effective 
velocity becomes appreciable compared to the speed of 
light. As a striking consequence, photons are radiated 
at frequencies higher than the mechanical frequency Q. 
A process of frequency up-conversion thus exists in the 
opto-mechanical coupling between vacuum fluctuations 
and mechanical motion of scatterers. 




FIG. 3. Spectrum of the radiation emitted by the cavity 
for Uos = 0.9 and a reflection coeflicient of r = 0.99. The 
top line of the frame corresponds to Qgjf/4. The peaks cor- 
respond to cavity resonance frequencies. The spectrum is 
plotted for a cavity oscillating globally at a mechanical fre- 
quency fl = 3n/L. The dashed line constitutes the envelope 
of the spectrum. It was obtained by averaging the spectrum 
of a single mirror over the effective velocities corresponding 
to different number of roundtrips. Photons are created at 
frequencies higher than the mechanical oscillation frequency 
through frequency up-conversion in the opto-mechanical cou- 
pling between vacuum fluctuations and the mirrors motion. 
Furthermore the radiation spectrum vanishes for frequencies 
equal to a multiple integer of the mechanical excitation fre- 
quency. 

A corresponding situation is found for a single oscillating 
mirror (dashed line) which is however not realistic as it 
would imply a mirror's mechanical velocity appreciable 
compared to the speed of light. The use of a cavity allows 
to reproduce the same spectral density within the band- 
width of the cavity resonance lines for realistic mirrors' 
velocities. A second striking feature is that no photons 
are emitted at frequencies equal to multiple integers of 
the excitation frequency f2 neither by a single oscillating 
mirror nor by a vibrating cavity. 

When comparing the cavity radiation spectrum to ex- 
pression ( ^tI) corresponding to a single mirror, a differ- 
ence is the emergence of peaks typical of cavity reso- 
nances. In fact, the interferences between the pathes 
corresponding to different numbers n of round-trips are 
essentially determined by the factors r'^^^'^™'^^'^ and 
j,n^-~2tT:Ku(n+i) ^ rpj^^ propagation dephasing after one 
round-trip is e'^'^^^'^ where K is the order of the mechan- 
ical frequency as compared with the fundamental reso- 
nance frequency of the cavity. It follows that the peaks 
are apparent at frequencies equal to an integer multi- 
ple of K~^, as shown on Fig. ^ with K = 3. Their 
shape is Lorentzian for a high finesse cavity. The width 
of each peak is given by the inverse of the cavity finesse. 
The number of peaks fitting into the interval [0, il] cor- 
responds to the order K of the excited cavity mode com- 
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pared to the mechanical frequency. 

VI. DISCUSSION 

In this paper we have presented a non-hnear calcu- 
lation of motion-induced radiation from a cavity taking 
fully into account the accumulation of dephasing through 
successive reflections of the field on partly transmitting 
mirrors. This approach confirms the main results of 
the linearized treatment which was previously used and 
makes it possible to specify its range of validity. Further- 
more the non-linearity due to the accumulative field de- 
phasing produces particular signatures of motion-induced 
radiation which cannot be calculated within the linear 
approximation. 

In the experimentally relevant case where the mirrors 
move with a velocity small compared to the speed of light 
the emitted photon number from a single mirror moving 
in vacuum scales with its squared velocity. Compared to 
this situation motion-induced radiation from an oscillat- 
ing cavity is enhanced by the cavity finesse. For high fi- 
nesse cavities as they exist for instance in the microwave- 
regime this enhancement brings motion-induced radia- 
tion within reach of an experimental observation. This 
clearly proves the cavity to be a much more favorable 
system for the generation of motion-induced radiation. 

In addition the present calculation shows that the lin- 
ear approach is valid when the effective rapidity, given 
by the mirrors' physical velocity multiplied by the cav- 
ity finesse, is much smaller than 1. We have given here 
expressions having a much larger range of validity. 

In order to measure motion-induced radiation it is 
necessary to dispose of signatures which permit to dis- 
tinguish vacuum radiation from spurious effects. The 
present calculation has allowed to identify two quantities 
showing signatures which could serve to this aim, the 
temporal variation of the radiated energy density and 
the spectral density of the emitted photon number. 

We have studied the emitted energy density as a func- 
tion of different effective rapidities. With increasing ef- 
fective rapidity the energy starts to concentrate in pulses 
which are emitted periodically into vacuum by the cavity. 
These pulses become the higher and the sharper the more 
the effective rapidity approaches its threshold value. The 
energy density diverges when the single-reflection rapid- 
ity equals half of the cavity losses during one roundtrip 
(a = p/2). The characteristic temporal variation which 
allows high energy densities in regularly spaced and nar- 
row time windows might be exploited in an experimental 
observation. 

The spectrum of motion-induced radiation shows sev- 
eral remarkable features. First photons may be radiated 
at frequencies higher than the mechanical frequency 
in contrast to the prediction of the linear treatment. 
A process of frequency up-conversion thus takes place 
in the opto-mechanical coupling between vacuum fluc- 



tuations and mechanical motion of scatterers. Second 
the spectrum always vanishes for all multiple integers of 
the mechanical oscillation frequency. Due to the opto- 
mechanical resonance condition motion-induced radia- 
tion is furthermore only predicted at particular frequen- 
cies corresponding to fractions of the mechanical oscilla- 
tion frequency. These signatures are different from pick- 
up effects and could serve to identify motion-induced ra- 
diation. 

So far we have discussed the behavior of the system in 
a regime where the cavity amplifies the dissipative effects 
of vacuum fluctuations. However, as a consequence of the 
divergence of the energy density there exists a threshold 
above which the system will show self-sustained oscilla- 
tions in analogy with an optical parametric oscillator [p^ . 
This regime is reached if the cavity finesse is increased 
above its threshold value. The amplification of motion- 
induced radiation should then exceed the cavity losses. 
It is to be expected that in this regime the cavity will 
emit photon pulses with much larger intensity than be- 
low threshold. If it were possible to reach this regime ex- 
perimentally an observation of motion-induced radiation 
as well as of its characteristics could be achieved more 
easily. It might thus also be interesting to calculate the 
radiated energy and the spectrum above threshold. The 
question then arises which are the mechanisms limiting 
the amplification of radiation in this regime. 

In conclusion, these results confirm the idea that it 
might be possible to show experimental evidence of the 
dissipative effects of motion in quantum vacuum. 
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